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Abstract 

The resuhs of measurements of the production of neutron-rich nuclei by the fragmentation of 
a ^'^Ca beam at 142 MeV/u are presented. Evidence was found for the production of a new 
isotope that is the most neutron-rich silicon nuclide, ^"^Si, in a net neutron pick-up process. A 
simple systematic framework was found to describe the production cross sections based on thermal 
evaporation from excited prefragments that allows extrapolation to other weak reaction products. 
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I. INTRODUCTION 



The study of properties of the most exotic isotopes continues to be one of the impor- 
tant tasks in experimental nuclear physics. In addition, masses, lifetimes, and properties 
of excited states are important not only for models of nuclear structure but also for the 
understanding of astrophysical processes. The first step in the study of a new exotic nucleus 
is its observation, which for neutron-rich nuclei demonstrates its stability with respect to 
particle emission. 

The neutron dripline is only confirmed up to Z = 8 (^^Oig) by work at projectile frag- 
mentation facilities in the US 1|, France , and Japan As indicated in Fig. [H the 
dripline rapidly shifts to higher neutron numbers at Z = 9 and ^^F22 has been observed by 
several groups . This shift makes the search for the neutron dripline in this region 

especially difficult but none the less important. Experiments at RIKEN in Japan [5] and at 



GANIL in France jol observed the two heaviest isotopes along the A = 3Z + 4 line, "'^Ne24 
and ^''Na26, by the fragmentation of ^^Ca28 projectiles. The heavier nuclei in this series, 
^°Mg28 and ^^AI^q, are unobserved at present. All nuclei with A = 3Z + 3 up to Z = 12 have 
been shown to be unbound. The neighboring nuclei with A = 3Z + 2 have been observed 
up to ^^Al28 but the production of the heavier nuclei from a ^^Ca beam requires a reaction 
with a net neutron pick-up. 

In the present work a series of measurements was carried out to search for new neutron- 
rich isotopes in this region and to measure the cross sections for production of these isotopes. 
A particular candidate for study is ^^Siao, a nuclide that has two more neutrons than the 
projectile nucleus. Nucleon pick-up products have been observed among fragmentation 
products, see for example Refs. js,!^], [lO], but their cross sections are significantly lower than 
those of pure fragmentation processes. The new data for the production cross sections builds 
upon the recent results from Mocko et al. and can provide a path to the production of 
the most neutron-rich nuclei. 



II. EXPERIMENTAL DETAILS 

A 142 MeV/u "^^Ca beam from the coupled cyclotron facility at the National Supercon- 
ducting Cyclotron Laboratory was used to irradiate either a ^Be target (724 mg/cm^) or a 
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FIG. 1: The region of the chart of nuchdes under investigation in this work. 

nat-^ target (1111 mg/cm^) located at the normal target position of the A1900 fragment sep- 
arator [ll|. The tungsten target was used due to its high melting point even though it is not 
monoisotopic. The average primary beam intensity for the measurements of the most exotic 
fragments was 70 pnA. The momentum acceptance of the separator was either Ap/p = ±1% 

ysis pro- 



or ±2% and the angular acceptance was 8 msr. The experimental setup and ana^ 
cedures used for this experiment were similar to those described in Refs. [3, H, [l2,ll3| and 
only the differences will be briefly described. The time of flight of each particle that reached 
the focal plane was measured in two ways: first, over the 17.8 m flight path between a plastic 
scintillator (22 mg/cm^ thick) located at the second dispersive image (image 2) and a 10 cm 
thick plastic backstop scintillator located at the focal plane of the separator, and also over 
the entire 35.6 m flight path of the A1900 fragment separator by measuring the arrival time 
relative to the phase of the cyclotron rf-signal. The magnetic rigidity for each particle was 
determined by the separator setting plus a correction based on the position measurements 
at image 2 with the plastic scintillator, and at the focal plane of the separator using a set 
of parallel-plate avalanche counters (PPACs). The standard focal plane detector setup was 
augmented to have three silicon PIN diodes (50 x 50 mm^ by 496 fim, 528 fim, and 526 fim 
thick) to enable multiple measurements of the energy-loss of the fragments and thus provide 
redundant determinations of the nuclear charge of each fragment. The simultaneous mea- 
surements of multiple AE signals, the magnetic rigidity, a scintillator signal proportional 
to the total energy, as well as the flight times for each particle provided an unambiguous 
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identification of tlie atomic number, cliarge state, and mass of the produced fragments. The 
position and angle measurements with PPACs at the focal plane also enabled discrimination 
against various scattered particles. 

The relative beam current was monitored continuously by a small BaF2 crystal mounted 
on a photomultiplier tube near the target position that provided a normalization for the data 
obtained at different magnetic rigidities. In order to map out the momentum distributions 
of the fragmentation products and provide the production yields, the magnetic rigidity of 
the separator was varied stepwise from 4.13 Tm to 4.89 Tm. The momentum distributions 
of isotopes between magnesium and phosphorus that were present at these settings were 
analyzed during the experiment. These measured distributions are in good agreement with 



liseH — h IJ] calculations, using either the Universial Parameterization 15j or the model 
by Morrissey so that the optimum separator setting for the heaviest isotopes (that 
were produced at a higher rigidity setting with very low rates) could be inferred from our 
LISEH — h calculations. Once the optimum setting was determined, an achromatic energy-loss 
degrader (^"^Al, 151 mg/cm^) was inserted at image 2 of the A1900 separator in addition to 
the plastic scintillator to cut down the range of atomic numbers of the fragments reaching 
the focal plane. 

A search for ^"^Si was carried out by performing several runs totaling 4.3 hours with 
the tungsten target and 5.8 hours with the beryllium target at a setting optimized for 
3^Mg and "^^Al at a rigidity of 5.045 Tm (4.9598 Tm after image 2). The combination 
of the higher energy loss of silicon isotopes in the thick targets and the image 2 degrader 
plus the expected momentum downshift due to nucleon pickup (cf. [9]) placed ^^Si in the 
acceptance of the fragment separator. Three events identified as ^'^Si nuclei were observed 
during the measurements with the tungsten target (see Fig. [2]) and none were observed 
with the beryllium target. The overall efficiency of the system was found to be 73^J5% and 
39l}2% when running with the tungsten and beryllium targets, respectively. The efficiency 
was dominated by the deadtime of the data acquisition system and the discrimination against 
pileup events in the focal plane detector. Trigger rates averaged 200 Hz for the runs with the 
tungsten target, and 450 Hz with the berryllium target. The simulated angular transmission 
ranged from 77% for '^'^Mg to 84% for ^'^Si with an estimated uncertainty of 5% using the 
technique described by Mocko et al. [lo| . 
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FIG. 2: (Color online) Particle identification plot of atomic number Z versus mass-to-charge ratio 
AjQ for Z = 7 to 15. 

III. RESULTS AND DISCUSSION 

The cross sections for the production of neutron-rich sihcon isotopes from this work are 
shown in Fig. [3] and given in Table [T] along with the cross sections recently reported by Mocko 
et al. for the reaction of "^^Ca with ^Be and ^^^Ta at the same bombarding energy. For 
the purpose of comparison we will consider the tantalum and tungsten targets as equivalent. 
The cross sections for reaction with the tungsten target are larger than those with beryllium 
by factors that range from approximately 2.5 at A = 38 to about 9 at A = 42, values that 
become significantly larger than the ratio of the geometric reaction cross sections 

a,{W) (AV3(Ca) + AV^(W))^ 
ar(Be) ~ (Ai/3(Ca) + yli/3(Be))^ 

The data show a smooth decline with increasing mass number (or neutron number) up to 
A = 42, and then a precipitous drop by about a factor of 110 for the two silicon isotopes with 
more neutrons than the projectile. The slope of the data compares well to the epax 2.15 
systematics [l3] although the data sets lie below the predictions. The epax parameterization 
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TABLE I: Cross sections for neutron-rich Mg and Si isotopes observed in this work. 



Isotope a (W target) a (Be target) 

(mb) (mb) 

36Mg (5±l)xlO-*^ (6j:^)xl0-'^ 

37Mg (9^3) X 10-8 (1.6l°:^)xl0-8 

38Mg (4±l)xl0-8 (4±l)xl0-9 

4iSi (1.3l°|)xl0-5 

42Si (9 ± 3) X 10-7 (9+^) X 10-8 

43Si (5 ± 2) X 10-9 (9l|)x 10-10 
44Si (7 ± 5) X 10-10 



describes the products of limiting fragmentation that occurs at high bombarding energies and 
only depends on the sizes of the target and projectile. Closer comparison of the prediction 
to the data shows that the cross sections for ^^Si from both targets are more suppressed 
than the average of the lighter isotopes. This is consistent with the idea that the most 
neutron-rich nuclei come from the decay of excited primary fragments that are themselves 
even more neutron-rich (and suppressed by the process of significant neutron transfer from 
the target at these bombarding energies due to momentum mismatch). 

Models of nuclear reactions used for counting rate estimates, like the intranuclear-cascade 
plus evaporation model [isj or abrasion-ablation in LISEH — h [19|] can not reproduce the low 
yields of the exotic nuclei observed in this study. Thus it is not possible to make reliable 
predictions for further work. As a starting point, the cross sections in peripheral two-body 
reactions have been analyzed in the framework of the Qgg systematics for a long time |20| . 
The central idea of the Qgg systematics is that the products are created in a statistical, 
thermal process and the cross section should follow the expression 

aiZ,A) = f{Z)exp{QjT) or In (a(Z, A)) cx Qgg , 

where Qgg is the simple difference between the mass excesses of the ground states of the 
product and reactant nuclei and T is an effective temperature that is fitted to the data. 
Such an ansatz is reasonable at low energies when the nuclei undergo a slow transfer process 
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FIG. 3: (Color online) The cross sections for production of neutron-rich silicon nuclei from Ref. [lo| 
and the present work. The data with yl < 43 (i.e., N < 28) are compared to the epax systematics 
for limiting fragmentation 

and for the observation of projectile residues from mass-asymmetric reactions where the 
bulk of the excitation energy is carried by the heavy partner. Over the years a number 
of measurements of light products at higher energies have found some agreement with this 
model as can be seen in Fig. H] (left panels) for the data from this study combined with the 
data from Mocko et al. [10]. The data for the most neutron-rich isotopes in each chain tend 
toward straight lines but the bulk of the data with the highest yields, highest precision, 
and lowest Q-values behaves very differently. It is important to realize that Qgg is most 
sensitive to the mass of the lighter fragment since the binding energy changes most rapidly 
with neutron and proton number in the low mass region. Previous studies that were used to 
develop the Qgg systematics relied on the analysis of the distributions of the light fragment 



from reactions in normal kinematics [20|. In the present work the lighter fragment is the 
target residue in the case of the beryllium target, whereas it is the projectile residue in 
the case of the tungsten target. The dominant factor in the exponential is then either the 
unobserved fragment (beryllium target, panel (b) of Fig. Hj) or the observed fragment (heavy 
target, panel (a) of Fig. Hj). 

Projectile fragmentation is usually not described as a two-body process, but rather as a 



sudden process that forms an excited prefragment followed by statistical decay. Charity 



23| 



has pointed out that the sequential evaporation of light particles from sufficiently excited 
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FIG. 4: (Color online) The variation of the cross sections for the production of neutron rich nuclei 
as a function of the two-body Q values (Qgg, left panels a, b) and as a function of the one-body 
Q value (Qg, right panels c, d), see text for details. Upper panels (a, c) show data for W (Ta), 
lower panels (b, d) for Be targets. Each symbol is labeled with the respective mass number. 
Data from the present work (below the dashed lines in each panel) were combined with data from 



Ref. [IQ:]. Solid symbols represent Q- value calculations 



open symbols based on the recommended values 
the most neutron-rich isotopes for each chain. 
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jased on the measured mass values, and 
. The lines represent exponential fits of 
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nuclei follows a general pattern that leads to a somewhat uniform distribution of final prod- 
ucts. This uniform distribution underlies the EFAX systematics. In the usual case neutrons 
are emitted preferentially from excited nuclei until the point at which the ratio of the widths 
for statistical emission of neutrons to charged particles, Fjv/r^, becomes small. Note that 
this expression includes neutrons and protons bound in clusters as described in Ref. 23 1. 
The individual emission widths, Tn and Tz, contain a number of factors but most of these 
factors approximately cancel in the ratio and the largest remaining term is an exponential 
of the difference between the neutron and proton separation energies, S^i and Spi 

Tn/Tz oc exp {Sp - Sn) ■ (1) 

The separation energies contain the masses of the daughter isotopes, thus, we can expect 
an exponential dependence of the yield on the mass difference between the daughter nuclei 
for proton and neutron emission in this model. The masses are not known experimentally 
for most of the very neutron-rich nuclei in this study. In an attempt to extract the average 
systematic behavior the cross sections are plotted as a function of 

Qg = ME{Z = 20, A = A8) - ME{Z, A) (2) 

in Fig. m (right panels), where ME{Z, A) is the mass excess in MeV. Qg is a function that 
compares the relative binding energies of all of the projectile fragments without regard to 
the target nucleus and is a plausible basis for comparison of products from a process that 
creates a small set of highly excited nuclei that then statistically populate all of the available 
mass surface. The figure shows that this function provides an excellent systematization 
of the data with each isotopic chain falling on a straight line. Moreover, the slopes or 
inverse temperatures decrease with atomic number and go from about 1.2 (Ar from Be) to 
a maximum of T ^ 2.5 MeV (Mg and Si from Be and Ta). The line from the production 
of magnesium isotopes can be extrapolated to predict a cross section of 0.04 ± 0.01 pb for 
"^^Mg, as yet unobserved. 



IV. SUMMARY 



The study of the production of the most neutron-rich silicon isotopes provided evidence 
for the existence of a new isotope, '^^Si, in a high energy reaction that requires the net trans- 
fer of two neutrons to the projectile. The decline of the cross sections for the production 
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of silicon isotopes with increasing mass number was found to parallel the EPAX parameter- 
ization but at a lower level, up to the point that neutron pickup intermediates begin to be 
important. The measured cross sections for nuclei with more neutrons than the projectile 
fall by approximately two orders of magnitude below a logarithmic extrapolation from the 
lighter isotopes. 

The variation of the cross sections for a large range of reaction products were considered 
in the framework of the well-known Qgg systematics developed for low-energy two-body 
reactions. Only the tails of the distributions had the expected linear dependence and the 
applicability of this model to projectile residues from reverse kinematical reactions is entirely 

questionable. On the other hand, all of the available data were shown to follow a very smooth 
systematic dependence, independent of the target, with the mass of the observed fragment. 
An extrapolation of the data using the new one-body Qg systematics indicates that a search 
for ^°Mg is feasible. 
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